Complex networks created by aggregation 
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We study aggregation as a mechanism for the creation of complex networks. In this evolution 
process vertices merge together, which increases a number of highly connected hubs. We study a 
range of complex network architectures produced by the aggregation. Fat-tailed (in particular, scale- 
free) distributions of connections are obtained both for networks with a finite number of vertices 
and growing networks. We observe a strong variation of a network structure with growing density 
of connections and find the phase transition of the condensation of edges. Finally, we demonstrate 
the importance of structural correlations in these networks. 
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I. INTRODUCTION 

Fat-tailed distributions of connections characterize the 
complex architectures of many real- world networks 0, |2, 
y, 13 ■ Several mechanisms may accounted for this form 
of degree distributions of networks. (Degree is the total 
number of connections of a vertex.) The most popu- 
lar concepts imply self-organization [ai H 13 • The self- 
organization mechanism is responsible for fat-tailed dis- 
tributions in a wide circle of evolving systems (see, e.g., 
Refs. 0,01), and not only in networks. 

Usually, a very particular preferential attachment ver- 
sion of the self-organization mechanism is discussed, so 
that highly connected vertices preferentially attract new 
connections [M B Q; but there are other possibilities. 
In this paper, we consider agglomeration as a compet- 
ing possibility. It is known that aggregation processes 
effectively generate power-law distributions (see, e.g., 
Ref. IlOl for an example). In networks, the analogue 
is the merging of vertices. By this mechanism, vertices 
accumulate their connections (agglomeration of edges). 
This increases a number of highly connected hubs and so 
gives a chance to arrive at a fat-tailed degree distribution. 

Evidently, the merging of vertices should take place in 
cellular networks (merging proteins) as well as in many 
other real- world networks. For example, in various net- 
works of economic relations, merging and splitting of en- 
terprises are basic elements of the evolution. The same 
is valid for networks of software components, electronic 
circuits, networks of relations between social groups, or- 
ganizations, institutions, and parties, networks of sub- 
jects, networks of notions, etc. Simple evolving networks 
with merging vertices have quite recently been simulated, 
Ref. [ij , and the generation of fat-tailed degree distribu- 
tion has been successfully demonstrated. (For a similar 
process in bipartite graphs, see Ref. T^l). 
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In the present paper we provide a comprehensive de- 
scription of the process of the creation of complex net- 
work architectures by the merge of vertices. It is impos- 
sible to obtain a uniform picture for all networks of this 
type. So, we describe a set of typical types of behav- 
iors by considering a line of basic network models, which 
can be studied analytically. These models may be gen- 
eralized in a natural way to include clustering and the 
condensation of clustering. Many other variations are 
also possible. 

All the models that we study in this paper generate 
fat-tailed degree distributions. We consider both non- 
equilibrium networks with a fixed number of vertices, and 
networks, where this number grows. We use the mean de- 
gree fc of a network as a relevant parameter. Then the 
variation of network architectures with k is essentially 
characterized by the "f{k) dependence. So, our main re- 
sults are presented in the form: exponent 7 vs. the mean 
degree k. 

In most of the networks in this paper, the evolution is 
due to two parallel processes: (i) the merging of vertices 
and (ii) random attachment of new vertices. However, 
we also discuss networks, where the second channel of 
the evolution is splitting (fragmentation) of vertices. The 
range of scenarios is wide, but in most of them we find a 
phase with the condensate of edges (in other words, gela- 
tion). Above some critical value kc of the mean degree, 
a finite fraction of edges is attached to a few vertices or 
to a single vertex. This condensation, unlike the situ- 
ation described in Ref. |T3|, occurs in the homogeneous 
networks. In the "normal phase", a degree distribution 
is of a power-law form, P{k) oc k~'* . Moreover, we ob- 
serve that, rather unexpectedly, even in the condensation 
phase, normal vertices have a scale-free degree distribu- 
tion. 

A resulting picture may be complicated by the presence 
of correlations, which is typical for non-equilibrium net- 
works. We demonstrate the importance of degree-degree 
correlations, in the most succinct of these network mod- 
els. The paper is organized as is follows. In Sec. ^ we 
describe the models and present in detail our results for 
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FIG. 1: Processes creating network 0. The labels r indicate 
that merging vertices are selected at random. "+" indicates 
an added vertex, g > 1 is a relative rate of the addition 
process. 



each of them. The complete final information can be ob- 
tained from this section. In Sec. lIIII we present details of 
our analytical calculations and simulations. 



II. MODELS AND RESULTS 
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FIG. 2; Processes creating network A. 
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FIG. 3: Processes creating network B. Here, m = 2. 



In this section we describe basic models of networks 
evolving due to aggregation processes and present our 
results. For sake of brevity we consider only undirected 
networks, i.e., networks with undirected edges. 



as well as its mean degree. Network becomes more and 
more dense with time. 

The resulting degree distribution is of a power-law (an 
asymptotics) form with exponent equal to 3/2: 



A. Network O 

This is the simplest model. The evolution starts from 
a given configuration of vertices and connecting edges. 
Loops of length one are allowed. At each time step (see 
Fig.HI: 

(1) 9^1 new bare vertices are added to the network. 

(2) Two randomly chosen vertices merge. 

Obviously, the final state of the network is a set of bare 
vertices plus a single vertex with all the edges (actually, 
loops of length one) attached. This is what we call the 
condensate of edges. 

B. Network A 

The (large) number of vertices of this network, N , does 
not change during the evolution. Initially, there is an 
arbitrary configuration of N vertices connected by some 
number of links. At each time step (see Fig. [2J): 

(1) A new vertex is added to the network. This vertex 
is attached to a randomly chosen vertex. 

(2) Two randomly chosen vertices merge. 

So, the result of the merging of vertices of degrees k' and 
k" is a vertex of degree k = k' + k" . The total degree 
of the network linearly grows, K{t) — K{t = 0) -I- 2f, 



P(fc) ~ fc-3/2 . 



(1) 



The main part of the distribution is stationary, but its 
low-degree part and a cut-off in the high-degree range 
change with time. This ensure the growth of the mean 
degree with time. Condensation of edges is absent. 

Note that we assume that this network is sparse. In 
principle, we allow multiple connections and 1-loops. 
However, we believe that in the sparse network regime, 
they are not important if we are not interested in the 
position of the cutoff of the degree distribution. 

This network was simulated earlier Jjj- O ur analytical 
results confirm the observations in Ref. IllL 



C. Network B 

This is a growing version of model A. At each time 
step (see Fig.|3l: 

(1) q > 1 new vertices are added to the network. Each 
of these vertices is attached to a randomly chosen 
vertex by m edges. 

(2) Two randomly chosen vertices merge. 

The total number of vertices now grows: N{t) = N{t = 
0) + {q~l)t. The total degree is i4:(i) = K{t = 0) + 2qmt. 
So, the average degree approaches the finite value k — 
2qm/iq~l). 

The condensation of edges is absent in this model. Un- 
like model A, the stationary degree distribution of this 
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FIG. 4: Processes creating network C. m = 2. 



growing network is a rapidly decreasing function. If, how- 
ever, the rate of the growth is low, q — 1 -^ 1, then the 
power-law dependence ^ is realized in the range of de- 
grees below a size- independent cutoff, k ^ iri/{q — 1)^. 



D. Network C 

In this growing network, at each time step (see Fig.^: 

(1) q > 1 new vertices are added to the network. Each 
of these vertices is attached to a randomly chosen 
vertex by m edges. 

(2) Simultaneously, a randomly chosen vertex merges 
with its randomly chosen neighbor, and the con- 
necting edge disappears. 

Rule (2) means the preferential choice: the second vertex 
in the network is chosen with probability proportional to 
its degree. The number of vertices and the total degree 
grow as N{t) = N{t = 0)+{q-l)t and K{t) = K{t = 0)+ 
2{qm— l)i, respectively. So, the mean degree approaches 
the value 



— „ qm — 1 

k9i 2- > 2m. 

q~l 



(2) 



It turns out that if the network is sufficiently dense, 
namely if 



k > kc = 2to(1 + ^yl - l/m) . 



(3) 



than a finite fraction of edges is condensed on a single 
vertex (or, maybe, on a few vertices). This takes place 
when the rate of the grows is low: 



q < qc = 1 + \/l - l/m . 
The fraction of edges in the condensate 
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FIG. 5: The fraction M of condensed edges in network G 
as a function of the k/{2m) (a) and of the growth rate q 
(b). Gurves 1, 2, and 3 correspond to m = 3, 6, and oo, 
respectively. 



the condensate in the limit of A: — > oo. (Note, however, 
that we consider a sparse network.) 

One can easily understand this condensation phe- 
nomenon. For the evolution of the number of edges in 
the condensate, Kh{t) = Mt, that is, the "macroscopic" 
number of edges attached to the hub, one can immedi- 
ately write the following equation: 



dKh 
dt 



K 



K-Kh 



{q-l)t 



-2 



(6) 



Here, K(t) is the total degree of the network, and the 
second factor on the right hand side of the equation is 
simply the mean degree fcn of "normal" vertices (i.e., the 
hub is excluded) minus 2. Indeed, according to rule (2) of 
the model, the probability that the hub will be chosen for 
the merging is Kh/K. Each act of merging, in average, 
increases the number of connections of the hub by /cn — 2, 
which explains the form of Eq. IJJ). Consequently, 



M = 



Akm + Am 2qra — q^m — 1 



k{k — 2m) 



qm — 1 



(5) 



behaves as M oc (k — kc) oc (qc — q) near the condensation 
point (see Fig.O. One can see that all the edges are in 



M 



M 



2{qm- 1) 



2{qm -1)- M 
q~l 



(7) 



One can see that this equation has a non-zero solution M 
[exactly of the form |SJ|] only if the growth rate parameter 




k/(2m), q 

FIG. 6: The 7 exponent of the degree distribution of network 
C as a function of k/{2m) and q. Curves 1, 2, and 3 corre- 
spond to m = 3, 6, and 00, respectively. Note that the depen- 
dences on k/{2m) and q are identical, but the ranges of the 
normal and the condensation phases are inversed [see expres- 
sions (|HJ and @]. The condensation takes place for k/(2m) > 
kc/{2m) = 1 -|- \/l — 1/m and q < qc = 1 + \J\ — 1/m. 



q is less than the critical value qc given by expression 
Q. Note that similar equations may be written for for 
networks with splitting vertices (see below). 

"Normal vertices" (i.e., with "microscopic" numbers of 
connections) have stationary power-law degree distribu- 
tions (asymptotics) both in the normal and in the con- 
densed phases. The 7 exponent of the degree distribution 
P{k) ~ fc-T is 



ikm 

7 = 2-h-= 



2 

k — Am 



(fc-2)(fc-2m) 



= 2- 



mq^ — 2(7771 + 1 
q{m — 1) 



>2 (8) 



in the phase without condensate [k < kc, i.e. q > qc) 
and 



Am — Akm 2qm — mq^ — 1 

= 2 + T-^ -r4 7T > 2 (9) 



2/c(m-l) 



{mq- l)(g- 1) 



in the condensation phase, where k > kc {q < qc)- Fig.lHI 
shows how the exponent of the degree distribution varies 
with the mean degree. Note that in both phases, near 
the condensation point, 7 — 2cx \k — kc\ oc I9 — 9c|. 

In the critical point, the degree distribution has the 
form: 



P(fc) 



1 



k'^ In k 



(10) 



Thus, this network is scale- free both is the normal and 
in the condensation phases. The power-law form of the 
degree distribution in the condensation phase is rather 
unexpected. 

Let us explain this remark in more detail. A close anal- 
ogy of the problem under consideration is the emergence 
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FIG. 7: Processes creating network D. m — 2. 



of the giant connected component in a growing network, 
where a phase transition with the Berezinskii-Kosterlitz- 
Thouless singularity takes place J4., J^ JJ., J/7j . In that 
case, the evolution equation for the size distribution of 
connected components is very similar to the evolution 
equation for the degree distribution in our case (see the 
next section). In this analogy, the giant connected com- 
ponent is analogous to the condensate of edges attached 
to a vertex, and the size distribution of finite connected 
components is analogous to the degree distribution in our 
case. The point is that the size distribution of connected 
components was found to be rapidly decreasing in the 
phase with the giant component (see Ref. [l5>|), while in 
contrast, in the present situation, the degree distribution 
is scale-free in the phase with the condensate. 

The paper [ij was mainly devoted to the simulation 
of the "static" version of quite similar network without 
multiple connections and loops of length 1, A^ = const 
with a growing number of connections, the sparse net- 
work regime. Scale-free degree distributions with expo- 
nents exceeding 2, without any condensation, were re- 
ported. We do not consider precisely this situation here, 
since we focus on stationary degree distributions. These 
take place in the growing network. 



E. Network D 

At each time step (see Fig. U)): 

(1) q > 1 new vertices are added to the network. Each 
of these vertices is attached to a randomly chosen 
vertex by m edges. 

(2) Simultaneously, the end vertices of a randomly cho- 
sen edge merge together, and this edge disappears. 

At first sight, this model is close to model C. The num- 
bers of vertices and connections grow in the same way, 
and the mean degree is the same, Eq. . The only differ- 
ence is the way in which the merging vertices are chosen. 
One can treat this merging process as transformation of 
random edges with their end vertices into single vertices. 

In fact, there is an essential difference since now the 
choice of vertices, in principle, depends on correlations 
between the degrees of the nearest neighbor vertices in 
the network. Let us compare models C and D once again: 



^- 9 





FIG. 8: The 7 exponent of the degree distribution of network 
D as a function of k/(2m) (a) and q (b). Curves 1, 2, and 3 
correspond to tti = 3, 6, and 00, respectively. The solid lines 
are obtained by neglecting degree-degree correlations. The 
dashed lines qualitatively show results of the simulation. In 
the condensation phase, exponent 7 does not change with k 
and q, and is close to 3. 



(i) The evolution of model C produces degree-degree cor- 
relations but is not governed directly by them, (ii) The 
evolution of model D depends on the degree-degree cor- 
relations and, in its turn, produces these correlations. 

Strict calculations taking into account degree-degree 
correlations should be rather cumbersome. Furthermore, 
in related networks, studied in Ref. [ll|, correlations have 
not been found. So, we applied a simplifying ansatz: we 
assumed that correlations may be neglected. 

It turns out that with this assumption, equations for 
the degree distribution have a reasonable solution only if 
the mean degree of the network is below some value: 



k <kc^ 2.204™ - 0.1115 + 0{l/m) 



(11) 



or, equivalently, q > Qc = 10.815 — 4.446/m. In this re- 
gion, our calculations provide a power-law degree distri- 
bution. The 7 exponent of this distribution approaches 
infinity at the minimal possible mean degree 2m (i.e.. 
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FIG. 9: The average number of vertices of degree k in simu- 
lated network D for g = 5, 7, 12 and m = 3. P{k) = Ji{k)/N. 
The values q = 5 and q = 7 imply the presence of the con- 
densate. The simulated networks contain up to lO** vertices. 




FIG. 10: The dependence of the mean degree of the nearest 
neighbors of a vertex on the degree k of this vertex, knn{k) 
in simulated network D for g = 5, 7, 12 and m = 3. The 
simulated networks contain up to 10* vertices. 



q — > 00) and near kc behaves as 



V^ 



(12) 



(see Fig. (S)). This is in sharp contrast to the behavior 
7(A:) near the critical point of network C (see Fig. O. 
Note the value 3 of exponent j ai k — kc- 

Above kc, with our assumption that the correlations 
are absent, the only solution was found to be patholog- 
ical. This proves the importance of the correlations in 
this network, which may be especially important in the 
condensation phase. 

For studying these degree-degree correlations, we re- 
sort to numerical simulations following the rules of model 
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FIG. 11: Splitting processes (examples): (a) splitting of a 
randomly chosen vertex into a pair of unconnected vertices; 
(b) splitting of a random vertex into a pair of connected ver- 
tices. In the simplest situation, each of possible resulting 
configurations occurs with equal probability. 



in such a way that all possible resulting configura- 
tions are realized with equal probabilities. 

(2) A randomly chosen vertex of degree k splits into a 
pair of vertices of degrees k' + k" = k + 2, connected 
by an edge. Again, we assume that all possible re- 
sulting configurations occur with equal probabili- 
ties. 

The number of these splittings per time step, in principle, 
may be not equal to 1. 



III. DERIVATIONS 



D. The degree distributions of the resulting networks are 
shown in Fig.|3 Power-law- like degree distributions were 
observed both in the condensation phase and in the nor- 
mal one. Fitting has given values of the 7 exponent 
slightly above 3 at the studied values of the parameter 
q. We obtained the dependence of the mean degree of 
the nearest neighbors of a vertex on the degree k of this 
vertex, k„n{k). Of particular note here is that A:„„(fc) 
has to be computed with care: the loops of length 1 are 
not to be considered for the vertex whose nearest neigh- 
bors are under study. However, we take into account the 
nearest neighbors with one loops. This may be especially 
important in the condensation phase. The obtained de- 
pendences knn{k) are shown in Fig. 1101 The main con- 
clusions are as follows: 

(1) the network is correlated, and the degree-degree 
correlations are strong; the correlations are of as- 
sortative type; one reason for this is the fact that 
for nodes with a high k the self-loops contribute 
strongly to the degree emphasizing such tendencies. 

(2) the correlations are present both in the phase with 
the condensate and without it. Fig. ^| demon- 
strates a non-monotonous dependence of fcrm(fe) on 
9; 

(3) the condensation phase transition takes place at 
higher values of the mean degree than fee, given 
by Eq. (P|l : 

(4) the observed values of the 7 exponent of the degree 
distribution are restricted from above. 



F. Splitting vertices 

Instead of the process of the random attachment of 
new vertices in models A-C, one can introduce another 
channel of evolution — splitting (fragmentation) of ver- 
tices. In this paper we only touch upon two possibilities 
(see Fig. ^| : 

(1) A randomly chosen vertex of degree k splits into a 
pair of unconnected vertices of degrees k' + k" = k 



In this section we describe details of our analytical 
calculations. We use an analytical technique similar to 
that for aggregation processes in more traditional sys- 
tems 18, 19] (for various aspects of the aggregation pro- 
cesses, see, e.g., Refs. OliflillE^). 



A. Network A 

The evolution equation for the average number 7V(fc, t) 
of vertices of degree k in the network A at time t has the 
following form: 



N{k, t + 1) = N{k, t) + 4,1 + ^N{k - 1, t) - liV(fc, t) 



1) ^ N{k',t)N{k'\t)-^N{k,t). (13) 



k'+k"=k 



N 



Here the average is over the statistical ensemble. (A ran- 
dom network is a statistical ensemble: a set of configu- 
rations with their statistical weights.) Three first terms 
on the right-hand side of Eq. l(T^ describe the process of 
the addition of a new vertex and the attachment it to a 
randomly chosen vertex. The two last terms describe the 
merging of a pair of randomly selected vertices. The fac- 
tor 2 in the last term is due to the fact that two vertices 
merge together. 

Note that we assume that our networks are large. So, a 
merging pair almost surely has no common nearest neigh- 
bors, and we can ignore the emergence of multiple con- 
nections during merging, if we do not interested in the 
cutoff region of the degree distribution. 

The total number of vertices in the network, N, is con- 
stant, so the evolution equation for the degree distribu- 
tion P{k, t) = N{k, t)/N has the form: 



dP(k,t) 
dt 



^Sk,i+P{k-l,t)-P{k,t) 



J2 Pik',t)Pik",t)-2Pik,t) 



(14) 



k'+k"=k 



(as is usual, asymptotically long times are considered). 
Assuming a stationary form of the degree distribution 
(except of a time-dependent cutoff and of, maybe, a low- 
degree part of the distribution) , we arrive at a stationary 
equation. The Z-transform of the degree distribution (a 
generating function) is 



i(z)^^z'=P(fc), 



(15) 



fc=0 



So, in a Z-transformcd form, we have 



N{k,t+l) = N{k,t)+qSk 
qm 



+ ^Nik-l,t) 



N{t) 



N{t) 



N{k,t) 



+ A^ E N{k',t)N{k",t)-^N{k,t). (20) 

'^ '' k'+k"=k 



N{ty 



Here, the number of vertices N{t) ^ [q — l)t. The 
evolution equation for the degree distribution P{k, t) = 
N{k,t)/[iq-l)t]is 



= n2(z) + (z-3)n(z) + z. 
The solution of this equation is 

1, 



(16) 



niz)^-[3-z^^{9-z){l-z)]. (17) 

This gives n(l) — 1, as it should be with J^k^i^) ~ ^■ 
This root of the equation is chosen, since it must be < 
n'(0) = P(l) < 1. The equation gives n'(0) = P(l) = 
1/3. 

Near z — 1, 



n{z) = analytical terms -f (1 



N3/2-1 



(18) 



This corresponds to the degree distribution P{k) ^ 
^-3/2^ Eq. ^, since the form of a Z-transform near z — 1 
and the asymptotics of its original are related to each 
other in the following way 



{q~\)[tdtP{k,t)+P{k,t)]=q5k,m + qmP{k-l,t) 
-qmP{k,t)+ Y^ P{k',t)P{k",t)-2P{k,t). (21) 



k' + k"=k 



Assuming a stationary degree distribution, in a Z- 
transformed form, this is 



= gz" + n^{z) -I- qmzn{z) - [q{m + 1) + l]n{z) (22) 

[compare with Eq. (|16() for q = 1]. The solution of 
Eq. H22|l is analytical at z = 1, so that the degree distri- 
bution is a rapidly decreasing function. Suppose, how- 
ever, that the rate q is close to 1. Then deviations from 
non-analytical behavior H18|) of n{z) are observed only in 
the range 1 — z < (q — l)^/m. This results in a size- 
independent cutoff fccut ^ m/(g — 1)^ of the power-law 
dependence P{k) ^ k^-^/^. 



Network C 



n{z ~ 1) = analytical terms + (1 
< — > P(fc>l)-fc-^, 



\7-l 



(19) 



if 7 is non-integer. The result is rather typical not 
only for aggregation processes but also for general non- 
equilibrium networks, where the mean degree linearly 
grows with time (see discussion in Ref. 01 )• 

Initially, we have assumed that the resulting degree dis- 
tribution is stationary. In principle, one can made a more 
general assumption, e.g., for brevity, P{k,t) = t'^f{k), 
where a is some exponent, and f{k) is an arbitrary func- 
tion of k. After the substitution of this form into Eq. 114() 
and Z-transformation, we find that the solution, n{z,t), 
depends only on z, that is stationary. 



B. Network B 

The evolution equation for the mean number of vertices 
of degree k in network B is 



The evolution equation for the mean number of vertices 
of degree k in this network is of the form: 



N{k,t+l) = N(k,t)+qSk 
qm 



+ ^N{k-l,t) 



Nity 



N{t) 



N{k,t) 



(23) 



^ V N{k\t)k"N{k'\t) 

N'{t)k ,,^kk 



k'+k" = k+2 



Nit) 



N{k,t) 



N(t)k 



N{k,t). 



Note that while deriving this equation, we did not assume 
the absence of correlations between the nearest neighbor 
vertices. Rule (2) of the model ensures that the second 
vertex for merging (a random nearest neighbor of a ran- 
dom vertex) is chosen with the probability proportional 
to the degree of this vertex. This is true irrespective of 
the presence of correlations in the network. 



For the stationary degree distribution, in a Z-trans- 
formed form, we have 



= qz"^ H =n{z)zn'{z) — =zn'{z) 

z'^k k 

+qmzn{z) — q{rn + l)n(z) . 



(24) 



[Note that the Z-transform of the convohition with 

j:k'+k"=k+2Pik')k"p{k")is 

^{n(z)zn' (z)-P(0)OP{0)~z(P(0)lP(l)+P{mP(0))} . 

z^ 

We have P(0) = 0, so that this expression simphfies.] 
Recall that 



/fc-2 
g= = > 1 

k — 2m 



(25) 



Equation 1)24(1 is the Abel equation of the second kind. 
In a canonical form it looks as 



[z^~n{z)]n'{z) = kq[z"'+'^ + {mz^-{m+l)z)n{z)] . (26) 

For the detailed analysis of a very similar equation for 
the size distribution of finite connected components in 
growing networks, see Ref. 1151] • Note, however, that in 
general terms, Eq. ((26|l has a different solution than that 
in Ref. [13. 

If P{k = 0) = 0, Eq. ^ implies 



n{z - 0) 



k{k~2) 



k (to + 1) - k{2m + 1) - 2m2 



(27) 



For obtaining a large k asymptotics of P{k), one must 
find the solution n{z) of Eq. H2()|l with the initial condi- 
tion (|27|l near z = 0. One can check that this solution 
arrives at 1 at z = 1: n(l) = 1. 

The linearization of Eq. (|26|l near z = 1 shows that 
n{z) linearly approaches z = 1, with a derivative n'(l), 
which satisfies the equation: 

h2 + n'(l)]n'(l) = M^[_2m + n'(l)] . (28) 
k — Zm 

The meaning of n'(l) is the mean degree of a vertex with 
degree k = o{N), that is the mean degree of a "normal" 
vertex. The solution of this square equation is 
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FIG. 12: The schematic view of the solution of Eq. II26|I 
in the normal (1) and condensation (2) phases. The dashed 
lines show the corresponding 1 — k{l — z) dependences. In 
the condensation phase, the slope of the solution at z = 1 is 
lower than k. 



where the expression for n'(l) changes its form, k^ 
Qc are given by Eqs. ||21l and Q, respectively. 



and 



At fc < fee, the slope n'{l) equals the mean degree of 
the network, k. That is, all the vertices of the network 
are "normal". Above kc (i.e., at g < qc) 'n'{l) = 2qm = 
2TO(fc — 2)/{k — 2m) < k. This means that a fraction 
M ~ [k — n'{l)]/k [see expression Q] is in the condensed 
state. The variation of the solution of Eq. H26|l with k is 
shown schematically in Fig. El 

The large-degree asymptotic behavior of the degree 
distribution of "normal vertices" is obtained by analyzing 
the form of the solution of Eq. H26|) near z = 1. We pass 
to new variables z — I — ^, n(z) = 1 — n'(l)^ + v{^) in 
Eq. H2t)|l . The inspection of the resulting equation shows 
that at k ^ kc, the contribution v{S,) is a power-law func- 
tion at small ^, w(^) oc ^", where exponent a is greater 
than 1. The substitution of this form into the equation 
allows us to obtain a. Using the correspondence (fT^ re- 
sults in the formulas © and (O for the 7 exponent of the 
degree distribution in the normal and condensed phases, 
respectively (see Fig. EJ. 

In the critical point, k — kc, the solution of the equa- 
tion for v{^) has a more complex form with an additional 
logarithmic factor (for more detail, see Ref. [laj)- The 
resulting solution n{z ^ 1) is 



n'(l) 



k — Am - 



4fcTO, — 4?n| 



2(k-2m) 



(29) 



That is, there is a point kc (and the corresponding qc), 
satisfying the equation 



n{z) = 1 — kc{l — z) + const 



1-z 



In [const (1 — z)] 



(31) 



k^ — AkcTn + Am — , 



(30) The asymptotics of the original of this Z-transform is 
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(32) 



Here the contour c is around 0, within the unit circle. 
The contour c is deformed to the contour c', which comes 
from — oo to +0 along the cut of In and then returns to 
— oo by the other shore of the cut. This deformation was 
made to ensure the decrease of the exponent. Thus, we 
have arrived at the degree distribution p()|l at the critical 
point. 

One can show that at /c ^ kc, in the large networks, the 
degree distribution has two regions, (i) The power-law 
dependence with exponents ^ and ^ is realized in the 

range of degrees In A; > 2 (A: — 4fc + 4m)/|fc —Ak'm + Am\. 
(ii) At smaller degrees, the critical dependence l(Tn|l is 
present. 



D. Network D 

Merging vertices in this model are the ends of randomly 
selected edges. In this case, for the strict description of 
the network evolution one has to solve an equation for a 
joint distribution P{k' , k") of the degrees of the nearest 
neighbor vertices. This is an essentially more hard prob- 
lern_than the analysis of Eq. (^H), |(5l]), and |(21l for P{k) 
or N(k) in previous sections. Instead of making these 
cumbersome calculations, we ignore the possibility of the 
degree-degree correlations and check whether this ansatz 
leads to reasonable results or not. This simplification al- 
lows us to consider a more simple evolution equation. 

Assuming the absence of correlations between degrees 
of the nearest neighbors of network D results in the fol- 
lowing rate equation: 



N{k,t+l) = N(k,t)+q5k^,n 

+ ^ _, Yl '^'^(^' ' ^) k"N{k" , t) 

N' {t)k k'+k" = k+2 



-2— ^7V(yfc,i). 
N{t)k ^ ' 



(33) 



For brevity, we consider only the case m > 2 (some of 
the resulting formulas may be different in the cases of 




FIG. 13: The solution of Eq.JSU_below (1) and above (2) the 
"critical" point fcc. Only for k < kc (i.e., q > qc), the solution 
provides n(l) = 1 and the limit slope k (see the dashed line). 



m — 1,2). So, for a stationary degree distribution, we 
have in a Z-transformed form: 



z-^k k 



+qmzn{z) — [q{ni + 1) — l]n(z) , 



(34) 



This is a nonlinear differential equation of the first order. 
It crucially differs from the corresponding Eq. (|24|l for 
the model C, and provides a different set of behaviors. 

Near z — 0, the solution of Eq. (|5Hl is of the following 
form: 



n{z - 0) 



k{k-2) 



2k{m - 1) - 4rn2 -|_ fc 



(35) 



This asymptotics is used as a boundary condition for 
Eq. ^. 

For a numerical analysis, the following form of Eq. H34|) 
is more convenient: 



n'{z) = kz-k{z^ - =-^ [k^{k - 2)2 



k — 2m 



+ {{k - 2)mz -k + 2m)n(z)] 



1/2 



(36) 



Note that only the solution with 'minus' is reasonable. 
One can check that at the boundary z = 1, if n(l) = 1, 
then Eq. l|M|) gives n'(l) = k and vice versa, if n'(l) = fc, 
then Eq. (|5Hl gives n(l) = 1. On the other hand, the nu- 
merical solution of Eq. H36|l shows that the solution with 
n(l) = 1 exists only at sufficiently low k. Above some 
value kc, the solutions of Eq. l|^ or Eq. l|^ turn out to 
be greater than 1 at some z: n(l) > 1 (see Fig. E| com- 
pare with Fig. [T^ . But this is impossible, since surely 
^^, P(fc) = 1. This contradiction indicates that our as- 
sumption does not hold at least for k > kc- 
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Let us study the analytical structure of the solution 
n{z), z near 1, for k < kc- As in Sec. 1111(71 introducing 
new variables, z = 1 — ^, n{z) = 1 — k^ + w(^), we pass to 
an equation for v{^). It turns out, however, that unlike 
Sec. IIII CI this equation has no solution v{^) ~ ^" with 
exponent 1 < a < 2. So, we have to search for the 
solution in the analytical form w(^) = C^^. Substituting 
this form into the equation gives 



C 



k{k + km — Am) 
4(fc - 2m) 



1-A 1-2 



(fc - 2){k - 2m)m{2k + to - 1) 
[k + km — 4to)2 



.(37) 



(One can check that only the sign "minus" at the square 
root leads to reasonable final values of the exponent of 
the degree distribution.) C is real only if A; < fcc, where 
kc is the point where the square root in Eq. (|37|l is zero. 
kc is exactly the point above which the solution n{z) 
is not reasonable [i.e., where n(l) > 1]. The resulting 
expression for kc is cumbersome but for to ^ 1 we have 



(38) 



which leads to the formula Hll|) . 

So, if fc < kc, we can substitute u(^) — C^^ + w{^) 
into the equation for v(^). An inspection of the resulting 
equation for w(^) shows that w{^) = D^^ , where expo- 
nent 6 > 2. One can easily find b: 




k (to 



1) 



(fc-2TO)(2C-fc) 



(39) 



with C given by Eq. H37|) . and by using the correspon- 
dence rule (|19|l we readily obtain the exponent of the 
degree distribution: 



7 



(fc - 2)(fc - 2m)m(2k -f to - 1) 



2W1-2 



(k + km — Am,)"^ 



fc(TO-l) / (fc-2)(fc-2TO)TO(2fc + TO-l) 



(40) 



fc+fc?7i — 4to 



(k + km^Am,)'^ 



E. Splitting vertices 

Here we only show the contributions due to splitting 
processes to evolution equations for the degree distribu- 
tion. Splitting process (1) of Sec. IIIFI generates with 
equal probabilities all possible configurations of two un- 
connected vertices [see Fig. ^J (a)]. For example, the 
splitting of a vertex of degree 4 produces: (i) with prob- 
ability 1/8, a pair of vertices of degrees and 4, (ii) with 
probability 4/8, a pair of vertices of degrees 1 and 3, and 
(iii) with probability 3/8, a pair of vertices of degrees 2 
and 2. One can easily write down the probability of re- 
sulting configurations for any degree of the splitting ver- 
tex. This allows us to find the probability that a vertex 
of degree fc will emerge due to the splitting of a vertex 
of degree q. Finally, instead, e.g., of the contribution 
Sk,i + P{k — 1, t) — P(fc, t) due to the process of random 
attachment of a vertex, we have the following terms: 



5:^(!)^('^'*)-^(^-'^) 



g=0 



(41) 



Note that actually only terms q > k are nonzero in this 
sum. In a Z-transformed form, this is 



2n 



1 + z 



t 1 — n(z, t) . 



(42) 



In splitting process (2) of Sec. IIIFI an emerging pair of 
vertices is interconnected by an extra edge. This leads, 
instead of the contribution of the form H41|l . to the fol- 
lowing terms: 



S2^(fc'>('^'^)-^(^'*)' 



q=t) 



(43) 



where only terms q> k ~ \ are nonzero in this sum. In 
a Z-transformcd form, Eq. (|43|l looks as 



2zn \ ——— ,t ) - n{z,t) 



(44) 



Formulas 14HI - (|44|I may be used to modify any of the 
evolution equations in the preceding sections. The result- 
ing equations, in a Z-transformed form, will be non-local 
due to the n{{l -\- z)/2) term. We do not analyze these 
functional equations in the present paper. 



This 7(fc) dependence is shown in Fig.|Sl 

The analytical results for network D were obtained 
in the framework of the simplifying ansatz: we ignored 
degree-degree correlations. The simulation of this net- 
work has shown that the correlations are significant (see 
Sec. IIIE|) and that our analytical predictions for the nor- 
mal must be corrected. Furthermore, the simulations 
have allowed us to describe the structure of the network 
in the condensation phase. 



IV. DISCUSSION AND SUMMARY 

Several remarks are necessary. We have studied a rep- 
resentative set of models allowing an analytical solution 
if no strong correlations are present. For sake of sim- 
plicity and brevity, we considered only models, generat- 
ing stationary (at least is some range of degrees) degree 
distributions. We have indicated that we had to ignore 
correlations in one of networks, and have demonstrated 
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that this neglect created some problems. In addition, 
we ignored multiple connections and loops of length 1, 
which, in principle, emerge during the evolution of these 
networks. 

The effect of multiple connections and loops of length 
1 in related network models was discussed in Ref. ^^ ■ 
Several features of the network structure depend on the 
presence (or absence) of these configurations of edges. In 
particular, the position of a size-dependent cutoff of the 
degree distribution may change. The number of vertices, 
attracting the condensate of edges may also change. We 
did not considered these details (cutoffs, the precise num- 
ber of vertices attracting the condensate, etc.). So, the 
strict accounting for multiple connections and one-loops 
should not essentially change the conclusions of this pa- 
per. 

One should emphasize the difference of the condensa- 
tion transitions considered in this work from those which 
were observed earlier in networks. We have shown that in 
networks C and D, scale-free degree distributions are re- 
alized both in the normal and in the condensation phases. 
In contrast, in equilibrium networks, a scale-free degree 
distribution exists only in the condensation phase, and 
a degree distribution is rapidly decreasing in the normal 
phase Hm. 

The merging of vertices is only one of processes in 
the networks under discussion. This process reduces the 
number of vertices, so other processes, increasing this 
number, must be present. In this paper, as a parallel 
process, we used an injection of new vertices (with sub- 
sequent attachment to existing vertices) and fragmenta- 
tion (splitting). This sub ject has been analyzed in the 
aggregation literature |lEl2a>l23- Also, one can imple- 
ment duplication (copying) or partial copying of vertices 
[la. 123 . l2g | or other processes. One should note, that 
correlations are important in low-dimensional aggrega- 
tion processes. In this respect we have an analogy with 
this case. In the usual aggregation theory the analy- 
sis proceeds by considering the scaling properties of the 
aggregation or fragmentation kernels (see e.g. |23| and 



references therein). For models such as those considered 
in our paper these are not known a priori and in the 
presence of degree correlations arise "self-consistently" . 
As is natural, a uniform picture for all networks cre- 
ated by aggregation processes cannot be obtained in prin- 
ciple. However, the models that we have considered in 
the present paper, reveal a number of basic features: 

(1) The aggregation easily generates network architec- 
tures where hubs play a profound role. 

(2) The aggregation often leads to gelation, or, in other 
words, to condensation of edges in these networks. 
We have found the condensation point at some 
mean degree value and have traced the variation 
of network structure in the normal and the con- 
densation phases. 

(3) These networks are evolving networks, and so their 
structure is characterized by strong correlations, in 
particular, by strong degree-degree correlations of 
assortative type. 

These features, which we observed by using demonstra- 
tive models, should be present in more complex networks 
of this kind. 
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